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ABSTRACT 

We study two nearby, early-type galaxies, NGC4342 and NGC4291, that host unusually massive 
black holes relative to their low stellar mass. The observed black hole-to-bulge mass ratios of NGC4342 
and NGC4291 are (!).9t.^%% and 1.9% ± 0.6%, respectively, which significantly exceed the typical ob- 
served ratio of ^0.2%. As a consequence of the exceedingly large black hole-to-bulge mass ratios, 
NGC4342 and NGC4291 are wS.lcr and ?a3.4cr outliers from the M, - Mbuigo scahng relation, re- 
spectively. In this paper, we explore the origin of the unusually high black hole-to-bulge mass ratio. 
Based on Chandra X-ray observations of the hot gas content of NGC4342 and NGC4291, we com- 
pute gravitating mass profiles, and conclude that both galaxies reside in massive dark matter halos, 
which extend well beyond the stellar light. The presence of dark matter halos around NGC4342 and 
NGC4291 and a deep optical image of the environment of NGC4342 indicate that tidal stripping, in 
which >90% of the stellar mass was lost, cannot explain the observed high black hole-to-bulge mass 
ratios. Therefore, we conclude that these galaxies formed with low stellar masses, implying that the 
bulge and black hole did not grow in tandem. We also find that the black hole mass correlates well 
with the properties of the dark matter halo, suggesting that dark matter halos may play a major role 
in regulating the growth of the supermassive black holes. 

Subject headings: galaxies: bulges — galaxies: elliptical and lenticular, cD — galaxies: individual 
(NGC4291, NGC4342) — galaxies: evolution — X-rays: galaxies — X-rays: ISM 



1. INTRODUCTION 

Observational studies of nearby galaxies have demon- 
strated a rather tight correlation between the mass 
of the central s upermassive blac k hole JM,) and the 
mass (Mbuige) (jMagorrian et al.l 119981 : iHaring fc Ria 



2004 ) or central velocity dispersio ii {a, 



20001 iFerrarese fc MerrittI l2000t iGiiltekin et al 



of the host spheroids. To explain the small intrin- 
sic scatter in the Af, — Mbuige and Af, — CTc rela- 
tions the theoretical paradigm has emerged, in which 
galaxy bulges and supermassive black hol es co-evolve 
and regulate each other's growth (e.g. 'Silk & Reesl 
1998; Fabian 1999; Kauffmann & Hachnch 2000; King 
20031 IDi Matteo et al.l l2005: Hopkins et al. 2006). How- 



ever, the masses of supermassive black holes also cor- 
relate well with the d ark matter halos of galaxies 
(jFerrarese fc MerrTtfl 120001 ) . suggesting that it may be 
the dark ma tter halo that (indire ctly) governs the black 
hole growth (IVolo nteri et al.l[2011| ) . Additionally, numer- 
ical studies indicate that the M, — Mbuigo relation can 
be reproduced by non-causal processes, which do n ot re- 
quire the co-evolution of bla ck holes and bulges ()Pens!J 
mM. [Tah^e fc Maccioll20TTh . 

In essence, despite the existence of scaling relations 
between the black hole mass and the physical properties 

|E-mail: abogdan@cfa.harvard.edu| 



of host bulges and halos, the origin of these relations 
remains a matter of debate. In this paper we employ 
X-ray, optical, and near-infrared observational data of 
two nearby early-type galaxies with unusually high black 
hole-to-bulge mass ratios to gain insight in the evolution 
of their bulges and black holes. 

The main f ocus of the work is NGC4 342, a lentic- 
ular galaxy (|de Vaucouleurs et al.l Il991[ ). which ex- 
hibits both an outer and a stellar nuclear disk 
(van den Bosch et al. 1998), and hosts a classica l bulge 
(ICretton fc van den BoschI 119991 iGraham et all 120011: 
iKormendv et al.l 120111) . The observed black hole-to- 
bulge mass ratio in NGC4342 is 6.9^3:3%, which signifi- 
cantly exceeds the typically observed '--^0.2%. Given the 
uncertainty of the black hole mass measurement and the 
scatter in the M, — Mbuige relation, NGC4342 is a wS.Ict 
outlier from the mean Af, — Mbuige relation (for details 
see Section 13. 3|) . Although the galaxy is projected on 
the Virgo cluster, it belongs to the W' cloud at a dis- 
tance of ~23 Mpc (Mci et al. 2007; Bogdan et al. 2012). 
Therefore we adopt a distance of D = 23 Mpc (1" — 111 
pc) for NGC4342. The Galactic colu mn density towards 
NGC4342 is iV^ = 1-6 X lO^" cm'^ (jDickev fc LockmanI 
[1993). 

To extend our stu dy we also investigate an ellipti- 
cal galaxy, NGC4291 (|de Vaucouleurs et aLlfTMl) . This 
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Fig. 1. — Left: 2MASS K-band image of a 6' X 6' (or 40 X 40 kpc) region around NGC4342. The center of the galaxy is marked with a 
white cross. The distribution of stellar light appears to be symmetric, thereby indicating the undisturbed nature of NGC4342. The main 
optical body is fairly small, ai90% of the stellar light is concentrated within a 0.7 X 2.2 kpc elliptic region. Right: 0.5 — 2 keV band merged 
Chandra image of the same region. One pixel corresponds to Ri2 arcsec. The X-ray emission has a very broad distribution, significantly 
exceeding that of the stellar light. The X-ray emission shows a distinct surface brightness edge and a tail, which we identify as a contact 
discontinuity or cold front — the nature of the edge is discussed in Bogdan et al. (2012). North is up and east is left. 



galaxy exhibits a black hole-to-bulge mass ratio of 
1.9% ± 0.6%, hence taking into account the black hole 
mass measurement uncertainty and the scatter in the 
M, — Mbuigc relation, the overly massive black hole cor- 
responds to a »3. 4o- o utlier from the mean M, — Mbuigo 
relation (Section 13. 3|) . NGC4291 is located in a poor 
group with 11 member galaxies (jGarcial 119931 ). For 
the distance of NGC4291 we adop ted D = 26.2 Mpc 
(1" = 126 pc) (jTonrvet al.l 120011) . The Galactic col- 
umn density towards NGC 4342 is Nh == 3.0 x 10^" cm"^ 
pickev fc Lockman|[T990l) . 

This paper is structured as follows. In Section 2, we 
introduce the data that we have used and describe its 
reduction. In Section 3, we compute the black hole-to- 
bulge mass ratios for NGC4342 and NGC4291. The ex- 
istence of massive dark matter halos around the galaxies 
is demonstrated in Section 4. Our results are discussed 
in Section 5, and we summarize in Section 6 

2. DATA REDUCTION 
2.1. Chandra 

NGC4342 was observed by the Chandra X-ray Ob- 
servatory in two observations (ObsIDs: 4687, 12955) 
for a total of 112.5 ks (Figure [l] right panel), whereas 
NGC4291 was observed in one pointing (ObsID: 11778) 
for 30.1 ks (Figure [5] right panel). The pointing for both 
galaxies is on the ACIS-S3 CCD, however for the analy- 
sis we used all available CCDs, namely ACIS-I2/I3 and 
ACIS-S1/S2/S3. 

Major ste ps of the data analys is were performed as 
described in lVikhlinin et all (|2005[ ). In particular, we re- 
processed the data set and applied the latest calibration 
corrections to the detected X-ray photons. We produced 
light curves for each observations to detect and remove 



high background periods (jMarkevitch et aLll2003[ ). After 
cleaning the flare contaminated time periods the clean 
exposure times were 78.8 ks and 26.2 ks for NGC4342 
and NGC4291, respectively. 

The diffuse emission fills the entire field-of-view of 
the ACIS-S3 CCD for both galaxies, therefore we de- 
termined the background level using "blank-sky" ob- 
servationfQ. Although the level of instrumental back- 
ground vari es slightly with time, its sp ectrum remains 
unchanged (jHickox fc Markevitchl[2006l ). To account for 
these variations we renormalized the background level 
using the 10 — 12 keV band count rate ratios, in which 
energy range the effective area of Chandra is virtually 
negligible. Since we aim to study the truly diffuse emis- 
sion from NGC4342 and NGC4291, emission from bright 
resolved sources need to be subtracted. Therefore, we ran 
the CIAO WAVDETECT tool in the 0.5 - 8 keV energy 
range to generate a source list. To obtain large source 
cells and exclude >98% of the source counts, we changed 
several para meters from their de fault values, which are 
discussed in iBogdan et al.l (|2012D . In total we removed 
126 and 43 sources from NGC4342 and NGC4291, re- 
spectively. The resulting source cells were used to mask 
out the location of bright point sources for further anal- 
ysis of the diffuse emission. 

2.2. Two-Micron All Sky Survey 

To study the stellar mass distribution we rely on the K- 
band images of the Two-Micron AU S ky Survey (2MASS) 
Large Galaxy Atlas (LGA) dJarrett et al.ll2003l ). which is 
known to be an excellent stellar mass tracer. The pub- 
hcly available K-band images of NGC4342 (Figure [T] left 

^ http://cxc.harvard.edu/contrib/maxim/acisbg/ 
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Fig. 2. — Le/i; 2MASS K-band image of a 5' X 5' (or 37.8 X 37.8 kpc) region around NGC4291. The center of the galaxy is marked with 
a white cross. The bright object to the east of the galaxy is a foreground star. Right: 0.5 — 2 keV band Chandra image of the same region. 
One pixel corresponds to ft!2 arcsec. The X-ray light has a much broader distribution than the K-band light. The diffuse X-ray emission 
appears to be more extended towards the south. North is up and east is left. 



panel) and NGC4291 (Figure [2] left panel) are not back- 
ground subtracted, therefore to estimate and subtract 
the background level we use several nearby regions off 
the galaxy. 

To convert the K-band luminosity to stellar mass, we 
compute the K-band mass-to-light ratios based on the 
B — V color in dices and r e sults of galaxy evolution mod- 
eling following HdleLaLl pOOl ): 

\og{M^/LK) = -0.206 -f (0.135 x {B - V)). 

For NGC4342 we use its luminosity weighted B — V = 
0.86 color index (Mihos et al. in preparation) that re- 
sults in Mi,/LK — .81, whereas in case o f NGC4291 we 
adopt B — V = 0.93 (|de Vaucouleurs et al.i 1991) yielding 
MjJLk = 0.83. The intrinsic scatter of the iBell et aLl 
(j2003| ) relation is ~0.1 dex, which determines the ac- 
curacy of the obtained mass-to-light ratios. The main 
sources of systematic uncertainties are the effects of dust 
and starbursts, which may generate ~25% uncertainty. 
The origin and signi ficance of uncert ainties is discussed 
in fuU particulars in IBell et all (|2003[ ). 

2.3. Deep optical image 

NGC4342 and its surroundings were observed as part 
of a deep, wide-field optical survey of the Virgo cluster 
by Mihos et al. (in preparation). The field surrounding 
M49 and NGC4365 was imaged in Spring 2011 using the 
Burrell Schmidt telescope at Kitt Peak National Obser- 
vatory. The imaging data consists of 73 1200 s images, 
each with a 1.65° x 1.65° field-of-view, dithered randomly 
by roughly half a degree between exposures. The indi- 
vidual images are then flattened using a master dark sky 
flat created from a median combine of 71 offset sky im- 
ages taken at the same hour angle and declination as the 
Virgo data. The extended wings of bright stars are mod- 
eled and subtracted from each image (.Slater ct al.,,2009, ). 



after which the images are sky subtracted by fitting a 
plane to the sky after masking out discrete sources. Fi- 
nally, the images are scaled to a common photometric ze- 
ropoint and median combined together to make the mas- 
ter mosaic with a native resolution of 1.45 arcsec/pixel. 
To bring out the faintest features, we mask individ- 
ual sources (stars and galaxies, both bright and faint) 
in the master mosaic, then rebin the image to 13 arc- 
sec/pixel using a spatial median filter. In this paper we 
use a 1 square degree cutout of this rebinned deep mo- 
saic, centered on NGC 4365. Because NGC 4365 lies 
near the edge of the Mihos et al. survey field, the total 
exposure time drops sharply across this portion of the 
mosaic, from 10 hours near NGC 4365 to 3 hours in the 
south-west corner of the field. More details on the imag- 
ing strategy and data reduction techniques for the V irgo 
imaging surv e y can be found in iMihos et al.l (|2005[) and 



imaging surv e y can 
IRudick et aH (^0H 



3. BLACK HOLE-TO-BULGE MASS RATIOS 

3.L M,/Mbuigo of NGC4342 

Stellar dynamical models of NGC4342 demonstrate 
that it harbors a supermassive black hole, whose mass 
(conv erted to 23 Mpc dis tance) is M, = 4:.6tll ^ 10** ^^o 
(iCretton &: van den Bosch 1999) . To estimat e the stellar 
mass o f the bulge we use the decomposition of lVika et alj 
(|2012| ). who found that ~ l/3 of the t otal s tellar light 
originates from the bulge. iVika et all ()2012[ ) computed 
that the apparent K-band magnitude of the bulge is 
10.3 mag, which corresponds to a K-band luminosity 
of Lk = 8.3 X 10^ Lk,0- To convert the luminosity 
to stellar mass we apply the K-band stellar-mass-to- 
hght ratio of Af^/i^ = 0.81 (Section [2?2|) . hence the 
bulge mass of NGC4342 is Mbuigc = 6.7 x 10^ M©. 
Thus, the black hole-to-bulge mass ratio of NGC4342 
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Fig. 3. — Black hole mass as a function of bulg e mass. Thick 
solid line shows the mean M, — A/buigc relation from lHaring fc Hixl 
II2004I) . whereas the thin dashed line represent the intrinsic scatter 
of the relation. Both NGC4342 and NGC4291 are highly significant 
outliers from the trend. 



is M./Mbuige = 6.9+^:1%. 

3.2. M./Mbuige ofNGC4291 

The supermassive black hole of NGC4 2 91 ha s been 
studied recently by iSchulze &: Gebhard^ (|2011[ ). who 
used stellar dynamical modeling to compute the black 
hole mass. They included the effect of a dark matter 
halo in their models, which resulted in a ma ss larger than 
previously reported (jGebhardt et al.ll2003() . Namely, the 
mass of the nuclear black hole in NGC4291 (for a dis- 
tance of 26.2 Mpc) is M , = (9.6 ± 3.0) x 10^ Mq 
(jSchulze fc Gebhardtl[20TTl ). The total apparent K-band 
magnitude of the galaxy is 8.417 mag, which translates to 
a K-band luminosity of Lk = 6.1 x 10^° Lk,©- Based on 
the K-band mass-to-light ratio of M^/Lfc = 0.83 (Sec- 
tion[221), the stellar mass of NGC4291 is Mbuigo = 5.1 x 
10^° Mq. Given these parameters the black hole-to-bulge 
mass ratio of NGC4291 is M./Mbuigc = 1-9% ± 0.6%. 

3.3. Comparison with the mean M, — Mbuigo relation 

The observed black hole-to-bulge mass ratios in 
NGC4342 and NGC4291 are unusually large compared 
to other galaxies. To illustrate this poi nt we sh o w the 
mean M, — Afbuigo relation from Hari ng fc Rixl (J2004D 
in Figure El with the loci of NGC4342 a nd NGC4291 
marked. Based on the iHaring fc Rixl ()2004D relation and 
the bulge masses of the two galaxies, the expected black 
hole masses are 7.7 x lO*' M© and 7.4 x lO'^ M© in 
NGC4342 and NGC4291, respectively. Thus, the ob- 
served values are factors of w60 and ~13 times larger 
than the predicted ones. From the intrinsic scatter of 
the relation (0.30 dex) and the uncertainty of the black 
hole mass measurements (0.18 dex and 0.12 dex), we con- 
clude that NGC4342 and NGC4291 are «5.1cr and «3.4ct 
outliers, respectively. 

Reversing the problem, we also compute the bulge 
masses, in which the supermassive black holes of 
NG C4342 and NGC4291 would be typical. According to 
the IHaring fc Rixl (|2004D relation and the observed black 
hole masses, we find that the black holes of NGC4342 and 
NGC4291 would be expected in bulges with 2.6 x lO" Mq 



and 5.0 x 10^^ M0 mass, respectively. These bulge masses 
exceed by factors of Ri39 and «10 the observed values in 
NGC4342 and NGC4291, respectively. 

4. TIDAL STRIPPING AS A POSSIBLE ORIGIN OF THE 
HIGH BLACK HOLE-TO-BULGE MASS RATIO 

One possibility to explain the unusually large black 
hole-to-bulge mass ratios observed in NGC4342 and 
NGC4291 is that most (>90%) of their stars were tidally 
stripped. However, the tidal stripping process would re- 
move not only the stellar content of galaxies, but also 
the more loosely bound dark matter halos. Therefore, if 
~90% of the stars were stripped from the galaxies, no 
significant dark matter halos should be observed around 
them. Thus, to test the stripping scenario it is critial 
to determine whether NGC4342 and NGC4291 host ex- 
tended dark matter halos. 

In Section 14.11 we use Chandra X-ray observations of 
the hot gas content of NGC4342 and NGC4291, to show 
that they reside in massive dark matter halos, thereby ex- 
cluding the stripping scenario. In Section [4. 2 1 we present 
a deep optical image of the surroundings of NGC4342, 
which - independently from X-ray observations - ex- 
cludes the possibility that majority of the NGC4342 stel- 
lar population was tidally stripped. 

4.1. Dark matter halos 

4.1.1. NGC4342 

The 0.5 - 2 keV band X-ray image of NGC4342 reveals 
a diffuse hot gas component associated with the galaxy, 
which exhibits a significantly broader distribution than 
the stellar light (Figure [T|). To compute the gravitating 
mass profile of NGC4342, we assume that the hot gas is 
in hydrostatic equilibrium (lMathewslll978HForman et aLJ 
119851: iHumphrev et al.ll2006l) and use the following equa- 
tion: 



Mtot(< r) 



kTgs^ (r)r I d In n^ d In Tgj 



G^m-p 



91ni 



9hn 



where Tgas and n^ are deprojected temperature and den- 
sity, respectively. To determine the projected tempera- 
ture and density profiles of the hot gas, we describe the 
soft-band emission with an optically-thin thermal plasma 
emission model (APEC model in Xspec). To obtain 
deprojected profiles, we use the technique described by 
[Ghurazov et al. (2003). Namely, we model the observed 
spectra as the linear combination of spectra in spherical 
shells plus the contribution of the outer layer. We assume 
that emissivity in the outer shell declines as a power law 
with radius at all energies. The matrix that describes 
the projection of the shells into annuli is inverted and 
the deprojected spectra are calculated by applying the 
inverted matrix to the observed spectra. 

Due the head-tail distributio n of the hot X-ray emit- 
ting gas (Figure [T] right panel: ' Bogdan et aLll2012D . the 
assumption of hydrostatic equilibrium is questionable at 
radii larger than '^5 kpc. To account for these uncertain- 
ties we computed mass profiles in three different sectors: 
(1) towards the northeast, (2) towards the southwest, 
and (3) assuming spherical symmetry. The left panel of 
Figure 2] illustrates that within the central 10 kpc re- 
gion of NGC4342 the hot gas is approximately isother- 
mal. Additionally, in .Bogdan et al., (,2012. ) we show that 
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Fig. 4. — Projected temperature profiles for NGC4342 (left panel) and for NGC4291 (right). For both galaxies the profiles are depicted 
for three different sectors: (1) sectors towards the north(east) are shown with blue, (2) sectors towards the south(west) are green, and (3) 
circular annuli assuming spherical symmetry are red. The corresponding position angles for the sectors are labeled in the figures. The 
position angle 0° corresponds to west and increases counter-clockwise. For NGC4342 the solid horizontal line shows kT = 0.54 keV, which 
was used to compute the gravitating mass profiles. To obtain the best-fit temperatures the abundances were fixed at 0.3 and 0.1 Solar for 
NGC4342 and NGC4291, respectively. 



the abundance is also fairly uniform within this region. 
Therefore to compute the mass profiles, we fix the gas 
tem perature at kT = 54 ke V, the abundance at 0.3 So- 
lar (jGrevesse fc Sau var'tOOS*) , and the column density at 
the Galactic value (Dickey & Lockman f990). We stress 
that beyond ~10 kpc the temperature of the hot gas 
is non-uniform, furthermore deviations from the hydro- 
static equilibrium significantly temper the accuracy of 
the mass measurements. Therefore, we rely only on the 
central 10 kpc region to estimate the gravitating mass 
of NGC4342. We note that within this region the total 
number of counts in the 0.5 — 2 keV band originating from 
the diffuse emission is 4713, which yields a signal-to-noise 
ratio of Ri50.8, allowing statistically accurate measure- 
ments of the mass profile. 

The mass profiles obtained for NGC4342 are shown in 
the left panel of Figure \5\ where we also show the la 
statistical uncertainties assuming that measurements in 
each radial bin are independent. In the same panel we 
also depict the stellar mass profile computed from the K- 
band luminosity using a mass-to-light ratio of M^,/ Lk = 
0.81, and the dynamical ma ss calculated assuming a 
circular speed of 220 km s^^ (jCretton fc van den BoschI 
Il999i) . Note that the circular speed was fixed in a cer- 
tain region (approximately 5" < r < 12"), where the 
contribution of random motions is not significant. The 
left panel of Figure [5] clearly demonstrates that beyond 
^1 kpc the stellar mass is significantly lower than the 
total gravitating mass. In particular, within the central 
10 kpc region the total gravitating mass is in the range 
(1.4 — 2.3) X 10^^ M0, which exceeds the stellar mass 
by an order of magnitude, implying the existence of a 



significant dark matter halo around the galaxy. 

Besides Chandra observations, the metallicity of the 
stellar population of NGC4342 indirectly and indepen- 
dently indicates that it formed in a massive dark matter 
halo. In the commonly accepted picture, galaxies with 
massive dark matter halos are able to retain a signifi- 
cant fraction of the a— elements produced by core col- 
lapse supernovae at early epochs of galaxy formation, 
whereas in galaxies with low mass hal os, a— elements 
are largely driven away in stellar winds (jTremonti et alJ 
120041: iGallazzi et all 120061 ) ■ At later epochs. Type la 
Supernovae enrich all galaxies with iron-peak elements, 
hence high a-to-iron metall icity ratios (or the oft en used 
[Mg/Fe] abundance ratios - iGonzalez et al.iri993D are in- 
dicative of massive dark matter hal os. NGC4342 bears a 
high [Mg/Fe] = 0.60 ± 0.04 ratio (jvan den Bosch et al.l 
1998), which is typical for galax ies with dark matter 
halos of a few times lO^^ Mq (jVazdekis et al.l 120041 : 
iHumphrev et al.ll2006( ). 

4.1.2. NGC4291 

In Figure [2] we show the K-band and the 0.5 — 2 keV 
band Chandra images of NGC4291. The truly difi'use X- 
ray emission associated with NGC4291 originates from 
hot X-ray gas and extends well beyond the optical light. 
To compute the gravitating mass profiles we assume that 
the hot X-ray gas is in hydrostatic equilibrium and follow 
the procedure described for NGC4342. Since the X-ray 
emission shows asymmetries, which may indicate devi- 
ations from hydrostatic equilibrium, we also computed 
the gravitating mass profiles using three different sec- 
tors: (1) towards the north, (2) towards the south, and 
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Fig. 5. — Gravitating mass profiles (Mtot) of NGC4342 (left panel) and NGC4291 (right) based on Chandra X-ray observations assuming 
that the hot gas is in hydrostatic equilibrium. We compute the gravitating mass in sectors towards the north(east) is blue, towards the 
south(west) of the galaxy is green, and assuming spherical symmetry is red. The solid lines represent the mean values, whereas the dotted 
lines indicate the lower and upper level of expected la uncertainties, which were calculated using bootstrapping method. The pink solid 
curves (Afstars) represent the stellar mass profiles obtained from the 2MASS K-band images. The black solid curve (M^y^) on the left 
panel shows an estimate of the total gravitating mass assuming flat circular speed of 220 km s~^ adopted from lCretton & van den BoschI 
JT999). Note that in the plotted regions the contribution of random motions is not significant. 



(3) assuming spherical symmetry. The temperature of 
the X-ray gas shows variations within the central 10 kpc 
region (Figure [3] right panel) , therefore when comput- 
ing the mass profiles, the temperature was left free to 
vary. The abundance of the hot gas is fairly uniform at 
0.1 Solar, therefore this parameter was fixed. Since the 
temperature cannot be determined accurately at radii 
larger than ~10 kpc, and the distribution of hot gas can- 
not be constrained at these larger radii, we only rely on 
the central 10 kpc region of NGC4291 to compute the 
gravitating mass profiles. Within the 10 kpc region we 
observed 3258 counts in the 0.5 — 2 keV band from the 
diffuse emission, hence the signal-to-noise ratio is high, 
«49.7. 

The resulting mass profiles with the la uncertainties 
are shown in the right panel of Figure [5l In the fig- 
ure we also plot the stellar mass distribution obtained 
from the 2MASS K-band image using the mass-to-light 
ratio of A/^/La- ~ 0.83. The plot illustrates that the 
stellar mass only plays an important role in the cen- 
tral ~2 kpc, beyond which radius the dark matter halo 
becomes dominant. Within 10 kpc the total mass is 
(2.1 — 2.6) X 10^^ M0, which is a factor of '--'5 times 
higher than the stellar mass of the galaxy. We thus con- 
clude that NGC4291 resides in a substantial dark matter 
halo. 

4.1.3. Excluding the tidal stripping scenario 

Although we cannot reliably characterize the dark mat- 
ter halo of NGC4342 and NGC4291 at radii much larger 
than ~10 kpc, Chandra measurements demonstrate that 



both galaxies host massive dark matter halos, which ex- 
tend well beyond the stellar light. 

In general, the stellar populations of galaxies are sig- 
nificantly more tightly bound than the associated dark 
matter halos. The dominant fraction of the stellar pop- 
ulation of galaxies is concentrated in the central ~5 kpc 
region, whereas t he dark matter halo rn ay extend up 
to a few 100 kpc (jHumnhrev et alllMffih . If ~90% of 
the stars were lost from NGC4342 or NGC4291, they 
were removed from the central few kpc region. However, 
owing to its low specific binding energy the dark mat- 
ter halo will be stripp ed first and faster than the stars 
(jLibeskind et al.l[2011[ ). Thus, the detection of a signifi- 
cant dark matter halo is inconsistent with the tidal strip- 
ping scenario. 

4.2. Tidal tail behind NGC4365 

The deep B-band image, shown in Figure [6l is also 
in conflict with the tidal stripping scenario. The image 
unveils an extended, faint tidal tail behind NGC4365. 
Interestingly, several galaxies, including NGC4342, ap- 
pear to be embedded in this tidal tail. NGC4365 lies 
0.33° (130 kpc in projection) northeast from NGC4342, 
and is presumably th e dominant member of the W 
cloud (jMei et al.ll2007f ). Although the exact origin of 
the tidal tail is uncertain, its presence and its distinct 
structure suggest that NGC4365 underwent a relatively 
recent (< 2 x 10^ years) interaction or merger with a 
companion galaxy. 

If -2 X 10" Mq stars were stripped from NGC4342 
during an interaction with NGC4365, a large fraction of 
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Fig. 6. — Portion of the deep B-band Virgo mosaic (Mihos et al. 
in prep.), showing a one square degree region around NGC4365. 
The mosaic has a native resolution of 1.45 arcsec/pix, but this 
image has been rebinned using a median filter to 13 arcsec/pix 
resolution to emphasize faint features. The limiting surface bright- 
ness in the image is ng ^ 29.5 mag/arcsec^, and regions brighter 
than ^g fa 26 mag/arcsec^ are rendered black. A prominent fea- 
ture on this image is a broad, ~200 kpc long tidal tail extend- 
ing southwest from NGC4365 with a mean surface brightness of 
/is Ri 28 mag/arcsec^ . A fainter and smaller tidal loop can also be 
seen to the northeast of NGC4365. 



the missing stars should be stih observable around the 
galaxy si nce NGC4342 h as an ~8 Gyrs old stellar popu- 
lation (v an den Bosch et al. 1998). However, the missing 
stars are unaccounted for, since the integrated B-band 
magnitude of the tidal tail behind NGC4365 (Figure [6]) 
is '^13.5 mag, which corresponds to Lb ~ 3.6 x 10^ Lb,0. 
We convert this luminosity to stellar mass using the B- 
band mass-to-light ratio of M^^/Lb — 3.6, which we es- 
timate under the (extreme) assumption that the tail has 
been stripped from NGC 4342, and using the B — V = 
0.86 color index of NGC4342 and results of galaxy evo- 
lution modeling (jBell et al.ll2003f) . Hence the total mass 
of the tidal tail behind NGC4365 is 1.3 x IQi" Mq more 
than an order of magnitude less than the "missing" stel- 
lar mass from NGC4342. Thus, the optical image, in- 
dependently from Chandra observations, disproves that 
~90% of the original stehar content of NGC4342 was 
stripped. 

5. DISCUSSION 

The exclusion of the stripping scenario implies that 
dominant fraction of the stellar population has not been 
lost as theorized in the tidal stripping scenario, but in 
fact never formed, making NGC4342 and NGC4291 gen- 
uine outliers from the Af, — Afbuigc relation. This conclu- 
sion has important consequences for the formation and 
evolution of black holes and galaxies. 

The presence of unusually massive black holes in two 
low stellar mass galaxies indicates that the growth of the 
supermassive black holes and the bulges of NGC4342 



and NGC4291 was not tightly coupled. In particular, 
our results imply that the growth of the supermassive 
black hole can precede the major burst of star formation. 
Although this concept is in conflict with the standard 
theoretical paradigm of bulge-black hole co-evolution, in 
which bulges and black holes co-evolve, a number of ob- 
servational and theoretical studies also support our con- 
clusions. In a recent study iReines et al.l ()2011f ) demon- 
strated the presence of a supermassive black hole in 
Henize 2-10 (a starburst galaxy without a substantial 
bulge) suggesting that supermassive black holes can form 
before their bulges. The evolution of the Af, — Afbuigo 
scaling relation at high redshifts also hint th at black holes 
may grow before their host bulges (e.g. iWalter et al.l 
[200l[Jahn"ke et al.|[2009l: iMerloni et al.ll2010[ ). Addition- 
ally, numerical simulations of quasar hosts indicate that, 
at high redshifts (5.5 S ^ ^ 6.5), black holes are more 
massive for their stellar hosts than exp ected from the 
relati ons observed in the local Unive rse (iKhandai et alj 
120111 ). In good agreement with this, iPetri et al.l ()2012f ) 
modeled the early formation of supermassive black holes 
and obtained similar results, namely that at high red- 
shifts the M, — Mbuigc relation was steeper than presently 
observed, hinting that black holes may grow faster than 
their host bulges. 

A major result of this paper is that both NGC4342 
and NGC4291 reside in massive dark matter halos. In 
fact, both the black hole masses and the observed dark 
matter halos are typical of galaxies having stellar masses 
that are ~10 — 40 times greater, hence the characteris- 
tics of the dark matter halos are consistent with those 
expected for the black holes. Therefore the only truly 
anomalous property of NGC4342 and NGC4291 are their 
low stellar masses. Since the black hole mass corre- 
lates well with the halo mass, it suggests that dark mat- 
ter halos may play a fundamental role in governing the 
black hole growth. Indeed, the tight relation between 
the black hole mass and th e dark ma t ter ha lo mass has 
already been re cognized by 'Ferr aresi (|2002D . Although 
iKormendv fc B ender (2011) hinted based on a sample of 
low- mass galaxies that supermas sive black holes do not 
correlate with dark matter halos, iVolonteri et al.l ()2011l ) 
re-examined their results and concluded that properties 
of dark matter halos do play an important role in the 
assembly of supermas sive black holes. The arguments of 
IVolonteri et al.l (|2011t ) are further supported by numer- 
ical studies, which demonstrate that the M, — Af bulge 
relation ca n be produced through t he central limit theo- 
rem (jPend 120071 : iJahnke fc Maccioll2011l) . hence the ex- 
istence of the scaling relation does not necessarily in- 
voke any causal process between the black hole and 
bulge evolution. Mo reover, numerical simulations of 
iBooth fc Schavd ()2010( ) pointed out that the primary pa- 
rameter in determining the black hole mass is the mass 
of the dark matter halo. 

Although NGC4342 and NGC4291 are outliers from 
the A/, — A'/buigc relation, both galaxies fit the M, — CTc 
relation rather well. The central velocity dispersions 
of NGC4342 and NGC4291 are 2 52.1 ± 8.4 km s~^ 
and 285.3 ±5.7 km s^^, respectively (jPrugniel fc SimienI 
119961 HyperLedsEI) . Based oi i the velocity disper sions 
and the mean M, — ac relation (|Giiltekin et al.l[2009[ ). the 

^ http://lcda.univ-lyonl.fr/ 
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expected black hole masses for NGC4342 and NGC4291 
are 4.2 x 10^ M© and 6.9 x 10* Mq, which are in good 
agreement with the observed values. This demonstrates 
that the velocity dispersion is a better pr oxy for the 
black hole mass than the bulge mass (see also lRusli et al.l 
120111) . Moreover, the tight relation between the black 
hole mass and the velocity dispersion is also in good 
agreement with the self-regul ated growth of supermas- 
sive black holes (ISi lk & Reel [1995 iFabianI [l999l: iKmal 
[200l IPTMatteo et al.ll20G5l ). 

The observed physical properties (M,/Mbuigo and 
M*/Mtotai ratios) of NGC4342 classify it as one of the 
most extreme outliers in the local Universe. Therefore, 
it is tempting to examine whether theoretical models 
can reproduce the observed characteristics of NGC4342. 
We study the galaxy c atalog constructe d on top of the 
Millennium simulation (jGuo et al.ll2010t ) and search for 
simulated counterparts of NGC4342. In particular, we 
searched for galaxies with M,/Mbuigo > 2.3% {2a lower 
than observed for NGC4342) and with bulge mass in the 
range of 2 x 10^ M© < Mbuigo < 1 x 10^° Mq. Several 
galaxies fulfilling these criteria are tidally stripped sys- 
tems, hence they do not host an extended dark matter 
halo. Among those galaxies with a notable dark matter 
halo, the highest observed black hole-to-bulge mass ra- 
tio is A/,/A/buigc ^ 3%, significantly lower than observed 
in NGC4342. Moreover, the simulated counterparts also 
host black holes with M, < 10* M©, and large fraction 
of the galaxies are strongly disk dominated. Thus, no 
galaxy in the Millennium simulation can reproduce the 
observed properties of NGC4342, thereby indicating that 
the evolution of NGC4342 is not fully explained by cur- 
rent theoretical models. 

Although the presently observed Af,/Mbuigc ratios of 
NGC4342 and NGC4291 significantly exceed those ob- 
served in other bulges, it is possible that the evolu- 
tion of these galaxies will eventually result more mas- 
sive bulges, hence normal M,/Mbuige ratios. Very low 
level star formation may be present in the central regions 
of NGC4342 and NGC4291 (Rafl^crty et al. 2008), which 
process could increase the bulge mass. However, low level 
star-formation (<IC 1 Mq yr~^) can produce the missing 
stellar mass, ^^(2.5 — 4.4) x 10^^ Mq, only in multiples 
of the Hubble time, moreover none of these galaxies host 
the required amount of either cold or hot gas to produce 
the missing stellar mass. Therefore, galaxy mergers may 
play a more important role in increasing the bulge mass 
of NGC4342 and NGC4291. However, the occurrence 
time scale of the sufficient number of (major-)mergers 
may be comparable or even longer than the Hubble time 
due to the relatively low merger rate a t the present epoch 
(e.g. iBowcr et al. 2006: iDe Lucia. G. fc Blaizot 20071 

Finally, based on the conclusions above we describe a 
potential evolutionary scenario, which is consistent with 
and could result in the observed physical properties of 
NGC4342 and NGC4291. The principal quantity de- 
termining the fate of the gas in the halo is the angular 
momentum. Indeed, large angular momentum promotes 
formation of the disk and more efficient star formation, 
but suppresses direct flow toward the central black hole. 
By contrast accretion onto a black hole and the coupling 
of the black hole feedback to the gas works best if the 
angular momentum of the gas is small (jChurazov et al.l 



[20051) . Thus, it is plausible that NGC4342 and NGC4291 
were formed in a sequence of accretion events, which 
brought mostly low angular momentum gas to the sys- 
tem. For details on the accretion of low angular momen- 
tum gas at the early-e pochs of galaxy formation {z ^ 6) 
refer to [Dubois et al.l ()2011.) . The accretion of low an- 
gular momentum gas resulted in enhanced black hole 
growth, while the star formation without strong rota- 
tional support mostly happened in a relatively compact 
central region. At intermediate redshifts {z ~ 2 — 4) 
the accreti on of gas with higher angular momentum 
([Dubois et al. 2011) could play a role in building the stel- 
lar body of NGC4342 and NGC4291. Once the black hole 
was large enough, it prevented further gas cooling (e.g. 
ISilk fc Rcci^'19 98HChuraz ov et al."2001^'Di Matteo et al.l 
2005; Hopkins etalj [2006; McNamara & Nulscn 2007), 
switching to the lo w accretion rate mode with very lit- 
tle star formation ([Churazov et al.l 120051 ). In this pic- 
ture, the virial temperature of the gas/halo is probably 
the most important parameter governing the black hole 
growth, hence it is not surprising that the final mass 
of the black hole "knows" about the characteristics of 
the halo. Therefore despite the significant deviations 
from the black hole-to-bulge mass relation, NGC4342 
and NGC4291 could still obey other relations, which are 
linked to their dark matter halo. 

6. SUMMARY 

We studied two early-type galaxies, NGC4342 and 
NGC4291, with black hole-to-bulge mass ratios of w6.9% 
and «1.9%, which are among the highest observed values 
in the local Universe. Our results can be summarized as 
follows: 

1. Assuming that the hot X-ray emitting gas content 
of the galaxies is in hydrostatic equilibrium we computed 
gravitating mass profiles. The total mass within 10 kpc 
is (1.4 - 2.3) X 10" Mq and (2.1 - 2.6) x 10" M© for 
NGC4342 and NGC4291, respectively. These values sig- 
nificantly exceed the stellar mass, implying the presence 
of significant dark matter halos around the galaxies. 

2. The detection of extended dark matter halos im- 
plies that tidal stripping, in which >90% of the stellar 
mass was lost, cannot explain the observed high black 
hole-to-bulge mass ratios. In case of NGC4342 a deep 
optical image of its surroundings further supports this 
conclusion. 

3. Since tidal stripping can be excluded as the origin 
of the high back hole-to-bulge mass ratios, we conclude 
NGC4342 and NGC4291 were formed as low stellar mass 
systems and that the black hole and bulge did not grow 
in tandem. In particular, our results suggest that the 
black holes suppressed star formation and the black hole 
possibly grew much faster than the host bulges. 

4. Since the characteristics of the dark matter halos 
are consistent with those expected from the black hole, 
we concluded that the dark matter halos can play a fun- 
damental role in governing the black hole growth. The 
black hole masses of NGC4342 and NGC4291 tightly cor- 
relate with the velocity dispersions, which indicates the 
self-regulated growth of the black holes. 

5. We described a potential evolutionary scenario that 
can produce the observed characteristics of NGC4342 
and NGC4291. In this picture, due to the accretion of 
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low angular momentum gas, the black hole grows more 
rapidly than the stellar population. Once the black hole 
reaches a critical mass, it prevents further gas cooling 
and star-formation. 
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